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Предложены оригинальные модификации генетических алгоритмов для расчета стационарных СОСТОЯНИЙ 
динамических систем химической природы, в том числе двиге ателей внутреннего сгорания, при заданных 
термодинамических параметрах. Особенностью разработки является разнесение малых и больших парамет- 
ров уравнений между уравнениями и алгоритмом их решения. Наряду с известными нейросетевыми методи- 
ками разработка интегрируется в состав имитационной системы «Мотор». 
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Введение. При изучении внутрицилиндровых процессов в поршневом двигателе внутреннего 
сгорания (ДВС) методами нелинейной динамики две практически важные задачи проблематично 
решить регулярными численными методами [1]. Первая состоит в отыскании стационарных и/или 
равновесных состояний, вторая — в быстром прогнозировании отклика эволюционирующей физи- 
ко-химической системы на управляющие воздействия. Источник трудностей — высокая размер- 
ность и многомасштабность рабочего процесса ДВС как динамической системы. Определенные 
надежды в предметных исследованиях возлагаются на генетические алгоритмы и нейросети [2]. 

Применительно к компьютерной имитации рабочего процесса в ДВС авторами разработа- 
ны оригинальные вычислительные методы: 

— генетический алгоритм отыскания равновесного и/или стационарного состояния физико- 
химической системы с реакциями и фазовыми переходами; 

— процедура оценивания решения задачи Коши применительно к впрыску, смесеобразова- 
нию, горению и образованию токсичных веществ в условиях ДВС. 
Постановка задачи. Определение равновесного состава некоторой фиксированной массы реа- 
гирующей, например, однофазной среды при заданных температуре и давлении, осуществляется 
путем минимизации суммарного термохимического потенциала системы [3] 


т, т; й 
б=У—^.| 9, +АТ пиру ти) —> тт (1) 
И; 7 
при материальных ограничениях 

от — т, =0. (2) 
В (1), (2) а -— массовая доля атомов го сорта в веществе } а = {о р ; т - вектор масс веществ, 


т = {ту}; т, - вектор масс атомов сорта, т, = {1}; и, И 9, - молярная масса го вещества 


а 
и его химический потенциал соответственно; Ми Р- общие масса и давление соответственно; 
9,(Г)=Н, + бр, (Т -298)-Т ($, + у Т / 298), где Н,, Ср,, $, - энтальпия образования, те- 
плоемкость при постоянном давлении и энтропия }-го вещества при температуре 7 = 298 К соот- 
ветственно. 





* Работа выполнена при поддержке Фонда АЕСОА. 
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Минимизированию (1), (2) посредством известных генетических алгоритмов [1] препятст- 
вуют жесткие ограничения (2): обычно вместо условного минимума (1), (2) ищут абсолютный ми- 
нимум модифицированного функционала 


И т) 2 : 
б=У—^ а, + КТ | п(МР) +т| —^ |||+(ат-т,) > тт. (3) 

ли, 7 
Метод решения. Классические генетические алгоритмы применяют для задач минимизации на 
компакте параметров: 


Е(х) > тт (4) 
аб ТЕМ, (5) 

Если в задаче имеются еще ограничения типа 
9(Х) <0 - покомпонентно, (6) 

их модифицируют следующим образом: 
9(х,о) = 9(Х)+5=0, с, >0, ]=1,2...М. (7) 
Задача сводится к 
Яка.) = од + о», [бя т (8) 
т 


с ограничением (5), где ^,>0, ]=12...М подбирается исследователем из расчета 


О] 


тит. Норму здесь берут евклидову |9,(х,5)|= (9,(х,5)}, а для подбора /, исполь- 
9,(Х, с, 








тах ([+(*] 


а,>х, >; 


пах (9,(%,с,)) 


а,>х,>В, 





зуют норму Чебышева, т. е. )., ^ 


При высокой размерности / ] задачи вероятность выхода популяции на ограничения (1) 
ничтожна, и процесс поиска приходится обрывать в одном из локальных минимумов. Понижение 
размерности системы путем исключения компонентов на основании (1) приводит к артефактам 
типа отрицательных концентраций или масс. 

Принципиальная новизна предлагаемого подхода состоит в переносе ограничений (2) из 
адаптивной функции в кодировку хромосома: в классическом алгоритме кодировка гена ограни- 
чена фиксированным интервалом [а, 6], а в модифицированном этот интервал приобретает опре- 
деленную подвижность в соответствии с ограничениями (2). При каждой генерации гена ищутся 
граничные значения аи В с учетом реализации предыдущего гена и ограничений (2). Физические 
соображения позволяют заметно сузить интервал локализации экстремума как сверху (ресурсом 
свободных атомов), так и снизу (устранением невязки материального баланса). Фактически алго- 
ритм позволяет выбрать базовые компоненты. Резюмируя, грубые локализующие свойства огра- 
ничений (2) выносятся в алгоритм, а тонкие уточняющие сохраняются в самих уравнениях. 

Если в классическом алгоритме кодировка гена ограничена фиксированным интервалом 
[а,6 | ‚ то в модифицированном этот интервал приобретает определенную подвижность в соот- 


ветствии с ограничениями (5), (6). При каждой генерации гена ищутся граничные значения а и 
Ь с учетом реализации предыдущего гена и действующих ограничений, т. е. для каждого го ге- 
на / =1.../\№М -1 ищется максимум при линейных ограничениях (10) и условии (11): 


т, > тах (9) 
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при линейных ограничениях 


о; + 9. >Т,.. +... +4, Ту = т), = Е, 


а, Т; + >, +... + ЧТ, = т, - Уа,,т, 
= 
7 (10) 


ОГП; + ЧТ, +. + ЧмуТь = Тм — У ант, 
= 
№М+1-/ 
и условиях 
1=1,2...М. (11) 


В результате находится верхний предел значения 17, , т. е. соответствующее значение РБ. 


т>т 


тм 7 


Генерация гена происходит также случайно, но уже с вновь вычисленным значением РЁ, т. е. 
Т, [а,6 (т, ) ‚ Где К =1.../-1. Соответственно первый элемент зависит только от ограничений 


(2). Решая на каждом шаге систему (9) — (11), например симплекс-методом, находим вектор т. 
Так формулируется универсальная методика нахождения минимума функционала с линейными 
ограничениями при помощи генетического алгоритма. Важно отметить, что система уравнений 
имеет высокую размерность только для первых т. 

Предложенная здесь модификация гарантирует выполнение различных сложных ограни- 
чений. При использовании в алгоритме гена большой длины (обычно применяют кодирование 
минимальной длиной 2 байта) и умеренного числа элементов этот алгоритм быстрее классическо- 
го. Максимальное число итераций для классического алгоритма Р„ и модифицированного Р, 


можно рассчитать по формулам 
Р, =2"М, (12) 
Р. =М?М, (13) 
где [ - длина гена, М - число ограничений, М - число элементов в системе. Из формул (12), 
(13) следует, что модифицированный алгоритм эффективен при № <2"?. При хорошей точ- 
ности М < 25/2 = 256 , а при высокой — М < 23? = 65536. 

Усовершенствованный алгоритм в сочетании со стандартным симплекс-методом решения 
(9)-(11) требует объемных вычислений при каждой генерации гена и таким образом оказывается 
неэффективным при решении задачи (1), (2). Поэтому здесь можно воспользоваться свойством 
системы а, > 0 и не решать (9)—(11), а просто ограничить следующий генерируемый ген остатка- 


ми масс элементов: 


11 
ти - аут, 
Ь, = мт 1 (14) 


1 
Ке[1,М] о; 





Фактически предлагаемый способ реализует выбор самим генетическим алгоритмом базо- 
вых компонентов в реагирующей системе, однако при этом не гарантируется строгое соблюдение 
материального баланса. Это несовершенство метода легко устраняется классической коррекцией, 
состоящей в замене минимизации (1), (2) на минимизацию (3). Таким способом локализующие 
масштабные свойства ограничений (2) вынесены в алгоритм, а тонкие уточняющие - в уравнения. 
В основу предложенного разделения масштабов положены естественные физические соображе- 
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ния: они позволяют заметно сузить интервал локализации экстремума как сверху — ресурсом сво- 
бодных атомов, так и снизу — устранением невязки материального баланса. 

При введении в систему свободных химических элементов (что не противоречит физиче- 

ским законам, но в некоторых случаях нежелательно ввиду их низкой концентрации) ограничение 
по минимальному остатку элемента удовлетворяет условиям (2) всегда: остаток элемента высту- 
пает как компонент системы. Такой подход является более общим при рассмотрении систем с 
большим числом элементов, компонентов и реакций. 
Реализация и результаты вычислительного эксперимента. Интерфейс компьютерной реа- 
лизации алгоритма показан на рис. 1. Термохимические константы реагентов, продуктов и про- 
межуточных веществ хранятся в проблемно-ориентированной БД. На данный момент также реа- 
лизованы функции импорта первичной термодинамической информации из стандарта программы 
$ТАМЗАМ [4]. 











[9 Генетический алгоритм = [5] х| 
ГНачальные условия ] Водород (Н2); 0.000858899824521248 
Вода (Н20 }; 0.0102699015793088 
Температура ТК) [222  — Количество итераций алгоритма Угарный газ (СО; 0.0159754024567025 
ПЕеРЫЯ 5 о 50 Эглекислый газ [С02]; 0.0188957961394675 
Исходная масса: 0.046 кг 
№ | [Название |Содинение| НН |5 [Ср [М | Сумма полученных масс: 0.046 кг 
| Водород Н2 0 130.6 28.83 0.002 [0.002 Массвобиранвиа на 00% 
|5 | М Вода Н20 -242000 188 30 (0.018 6 _ 
|8 | м НИЕ ыы Равновесный состав 





Эглекислый газ (С02 0.046 


Пуск алгоритма Графики 


а 
|= 
Я 












































0.04 


Исходный состав 


б00-Ч + 


Рис. 1. Результат расчета равновесного состава в системе 
СО)? (1 моль) + Н? (1 моль) <> СО + Н?2О при Т = 2222 К, Р= 7,5 атм 


Верификация алгоритма осуществлялась на интересных для авторов системах: крекинг уг- 
леводородов, горение традиционных нефтегазовых и перспективных спиртовых моторных топлив, 
продукты сгорания различных ракетных топлив. Пример такого сопоставительного расчета при- 
веден на рис. 2: посредством усовершенствованного генетического алгоритма определен равно- 
весный состав системы Н2О - Н> — ОН - О> — О -Н при давлении 5 атм в диапазоне температур 
2000-6000 К. Результаты сопоставления рассчитанного равновесного состава продуктов сгорания 
различных моторных и ракетных топлив с данными работ [5], [6] свидетельствуют, что главные 
компоненты химического равновесия определяются предложенным методом вполне надежно. 

Во всех изученных системах содержание основных компонентов определялось при помощи 
усовершенствованного генетического алгоритма с приемлемой погрешностью. Зная содержание 
всех ведущих компонентов, нетрудно рассчитать также важные, например, в плане токсичности, 
малые составляющие равновесной смеси. Это ценно, поскольку размерность задачи существенно 
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Доля 
вещества 


0.6 зависит от того, насколько верно удается априорно выбрать но- 
менклатуру главных компонентов в химическом равновесии. 
Перспективы развития биоинспирированных методов. 
Отрабатываемый авторами нейросетевой алгоритм быстрого 
оценивания решения задачи Коши сейчас «обучается» на фраг- 
ментных задачах: 

— имитация динамики распыла — испарения в условиях 
дизеля; 

— прогнозирование мощности и экономичности рабочего 
цикла по входным параметрам; 
он — определение содержания токсичных веществ в отрабо- 

Ри `  тавших газах. 
Рис. 2. Рассчитанные парциальные Следует отметить, что если «обучение» возможностям 
концентрации компонентов: первой и второй задачи удается ограничить расчетными «обу- 
ЗА НЕ Ньв- ОНУ 05 чающими» данными, для уверенного прогнозирования токсично- 
п — Н, У - О (сплошные линии — ы 
ани даннне Г) сти выхлопа приходится привлекать экспериментальный мате- 
риал. Тем не менее, об адекватности этого алгоритма целям со- 
вершенствования программного комплекса «Мотор» [7], [8] свидетельствует успешность «обуче- 
ния» нейросети на результатах представительного множества численных расчетов. Сеть позволя- 
ет прогнозировать токсичность выхлопа известного ей двигателя в пределах обучающей выборки. 
Однако в силу существенной нелинейности внутрицилиндровой физико-химической кинетики на- 
дежность такого прогнозирования пока недостаточна. 

Внедрение новых биоинспирированных имитационно-аналитических средств в состав сис- 
темы «Мотор», интерфейс и некоторые возможности которой показаны на рис. 3, расширяет ее 
функциональность, прежде всего, в части приспособленности для комплексной оптимизации и 
управления внутрицилиндровыми процессами. 
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КОНСТРУКТИВНЫЕ ПАРАМЕТРЫ ДВИГАТЕЛЯ: 











ри Число . [2] Объем ‚,_ [^ [^ 
Тактность {= ЕН цилиндров -и Е кс \“ 17000 5] ДП 2-8 
Щ— РЗ] Степень -® В 0.25 [м 
МаСТОТа| =- (9% сжатия а ПИ РЕКИ |=— 0.25 

















Г 
В А 
ид топлива гр. 6е+41 Осв,Дж } ЗОсь Ват 
Бензин Е 2е+71 
5е+4 
Состав топлива. С Н 1.5е+7 
‚. 4е+4 ь 
массовых % 
3е+4 
ПАРАМЕТРЫ РАЕ ный 


2е+4 
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-180 -120 -50 0 60 120 180 -180 -120 -50 0 60 120 180 
—— Тепловыделение в реакциях Темп тепловыделения в реакциях 











* Газодинамическая 
Расчетные энергетические характеристики двигателя 


Тепловыделение при 5 = 0 - 6.5% Р‚„„- 65.6 атм 
Тепловыделение при Р""- 95.8% Т,- 2657 К 
Длительность сгорания 9’ - $’ - 2.1 рад. ПКВ = 180 град. ПКВ ^ 
Мощность И’- 1624.9 л.с. КПД тепловой л-24.5 % Ту -180 град. пкв ^ 1388 К 
Удельный расход топлива - 245.3 п/л.с.ч. Я 14.2 атм 


- 115.7 Кград.ПКВ 


-Модель горения: 


С Химкинетическая << Вернуте 


4.9 атм 


пэх 


Жесткость сгорания се 


пох 


- 3.068 атмирад.ПКВ Ав 


Рис. 3. Интерфейс и функциональность программного комплекса «Мотор» [8] 
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Заключение. Предложено усовершенствование генетического алгоритма для расчета стацио- 
нарных состояний динамических систем химической природы, например, двигателя внутреннего 
сгорания. Особенностью разработки является разнесение малых и больших параметров уравне- 
ний между уравнениями и алгоритмом их решения. 

На данный момент интегрированная в комплекс «Мотор» [7], [8] система моделей и мето- 
дов, в том числе биоинспирированных, позволяет: 

— рассчитывать энергетические характеристики ДВС и токсичность выхлопа по СО, СН,, 
ВСНО, МО при использовании различных топлив и на переходных режимах; 

— выявлять ведущие физико-химические процессы во внутрицилиндровой кинетике; 

— анализировать динамику внутрицилиндровых процессов и вырабатывать способы их оп- 
тимизации; 

— формулировать критерии параметрической оптимизации; 

— разрабатывать эффективные способы управления. 
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АРРЫСАТТОМ ОР ВТО-ТМ$РТВЕО ТЕСНМТОЧЕ$ ТО МОБЕЕ ТОХТС РАКАМЕТЕК$ 
ОЕ 1МТЕВМАЕ СОМВИУ$ТТОМ ЕМСТМЕ 


0.5. СТВЗ$Н, Е.М. ГАДОЗНА, $.М№. КНОЕОООМА, 0.5. ТЗУМВАЕО\У 
(Роп Зае Теспткса! Упмег®Ку) 


$оте опдта!/ то!йсайоп$ оЁ {пе депейс а/догйт5 Юг апа!у2тд {те чаНопагу 55 о! ше спетка!/ дупатс 
зу${етс аЁ {Те аейпей {Пегтодупатис рагатеегз, ГС епдтез атопд о{ТЕГ$, аге опегед. Тре деуеортепе Геа- 
{иге 15 аШосаНоп оЁ зта!! апа !агде едиайоп рагатеег$ атопд {Те едиайоп$ апа пе! 5о/ийоп5. Аопд ий" пе 
Кпоит пеигопе! теосб, {Пе деуе/ортепе сап Бе педгайеа т 'Епдте’ ятшайоп зузет. 

Кеуигогаб&: депейс а/дог"т5, пеигопеЁ, сотризйоп, ехпаизЕ, юЮхсйу. 
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